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Abstract  {B^}  rotor  modal  displacement  of  (Y^,  y^^) 

A  newly  developed  global  dynamic  model  was  [C|j]  bearing  damping  matrix 

used  to  simulate  the  dynamics  of  a  gear  noise  rig 

at  the  NASA  Lewis  Research  Center.  Experi-  [CJ  [$] 

mental  results  from  the  test  rig  were  used  to  _ 

verify  the  analytical  model.  In  this  model,  the  [Cy  [^cl 

number  of  degrees  of  freedom  of  the  system  are 

reduced  by  transforming  the  system  equations  of  [CJ  casing  structure  damping  matrix 

motion  into  modal  coordinates.  The  vibration  of 

the  individual  gear-shaft  systems  is  coupled  [CJ  [CJ  J 

through  the  gear-mesh  forces.  A  three- 

dimensional  bearing  model  was  used  to  couple  the  {D}  rotor  modal  displacement  of  Z 
casing  structural  vibration  to  the  gear  rotor 

dynamics.  The  system  of  modal  equations  is  casing  modal  displacement  of  (Z^, 

solved  to  predict  the  resulting  vibration  at  sever¬ 
al  locations  on  the  test  rig.  Experimental  vibra-  {D^}  rotor  modal  displacement  of  0^ 

tion  data  were  measured  at  several  running 

speeds  and  were  compared  to  the  predictions  of  F  gear  force 

the  global  dynamic  model.  There  was  excellent 

agreement  between  the  analytical  and  experimen-  {F(t)}  external  and  mass  imbalance 
tal  vibration  results.  excitations 
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Nomenclature 

rotor  modal  displacement  of  (X, 


{F^(t)}  force  acting  on  casing  structure 

{FQ(t)}  nonlinear  gear  forces,  through  gear 
mesh  coupling 


casing  modal  displacement  of  (X^,  X^^)  {F,(t)}  shaft  bow  force  =  [KJ  {W,} 


rotor  modal  displacement  of  (Y,  0^)  [G^  rotor  angular  acceleration 
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[KJ 

casing  structure  stifbiess 

[KJ 

shaft  bow  stiffness 

m 

mass  matrix  of  rotor 

(MJ 

mass  matrix  of  casing  structure 

R 

pitch  radius  of  gear 

{[W]} 

generalised  displacement  vector  of 
rotor 

{[WJ} 

generalised  displacement  vector  of 
casing 

{W,} 

genearlised  displacement  vector  of 
rotor 

X,Y,Z 

displacement  vectors 

a 

angle  of  orientation 

e 

Eq.(5) 

0 

generalised  displacement 

torsional  rotational  vector 

lateral  rotational  vectors 

It  friction  coefficient  between  the  gear 

tooth  surface 

orthonormal  mode 

[$  J  orthonormal  mode  of  casing 

[it;]  critical  speeds  of  each  rotor 

[wj  critical  speeds  of  casing 

Introduction 

Large  vibrations  in  gear  transmission  systems 
causes  excessive  wear  and  crack  formation  in  gear 
teeth,  which  results  in  premature  gear  feulure. 
With  the  need  for  higher  operating  speeds  and 
power  from  transmission  systems,  the  problem  of 
excessive  vibration  becomes  even  more  critical. 

In  order  to  insure  smooth  and  safe  operation,  it  is 
necessary  to  understand  the  dynamics  of  the  gear 
transnoission  system. 

Two  areas  of  research  in  the  dynamics  of  gear 
transmission  systems  are  (1)  analytical  Simula* 
tion  and  (2)  experimental  testing.  There  is  a 
great  deal  of  literature  on  the  vibration  analysis 
of  a  single  gear  stage.^**  Some  work  has  been 
done  on  multistage  gear  vibration,^*^  but  very 
limited  work^*^^^  has  been  done  on  the  dynamic 
analysis  of  gearbox  vibration.  Considerable  effort 
has  been  devoted  to  experimentally  studying  gear 
dyanmics^^'^^  and  loc^ized  vibration  effects  on 
gear  teetb.^^  A  few  studies  have  been  conducted 
to  correlate  analytically  predicted  and  experi¬ 
mentally  measured  gearbox  vibrations. 

This  paper  correlates  the  experimental  results 
obtiuned  from  the  test  rig  at  the  NASA  Lewis 
Research  Center  with  predictions  from  an  analyt¬ 
ical  model  developed  by  using  the  modal  synthe- 
ris  method.^  The  major  excitations  of  the  rotor 
system  include  mass  imbalance,  shaft  residual 
bow,  nonlinear  gear-mesh  forces^*  and  gear-tooth 
frictional  effects.^^  The  vibratory  motion 
between  the  rotor  and  the  casing  is  coupled 
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through  the  support  bearing  in  the  iateral  and 
axial  directions.  ^  Gearbox  mode  shapes  and 
vibration  predictions  from  the  analytical  model 
are  compared  to  those  obtained  from  experimen¬ 
tal  testing. 


Analytical  Procedure 


Development  of  Equations  of  Motion 

The  equations  of  motion  for  each  gj^-shaft 
system  can  be  written  in  matrix  form'  as 

IM]  {w}  -H  IGJ  (w)  -h  [G^]  {W}  +  [Cfcl  {W  -  W,) 
+  [Kbl{w-w4  +  (K,]{w-  wj 
=  (F(t)}  +  (FcCt)} 
(1) 


where  [M]  is  the  mass  matrix  of  the  rotor  (iner¬ 
tia),  {w}  is  the  generalized  displacement  vector 
consisting  of  the  three  displacement  vectors 
X,Y,Z  with  the  corresponding  lateral  0j,  and 
torsional  0^  rotational  vectors  as 


{w}  = 


(X) 

(«x) 

(Y) 

(V 

(Z) 

(«t) 


(2) 


and  where  [G J  is  the  gyroscopic  force;  [GyJ  is  the 
rotor  angular  acceleration;  [Cj,]  is  the  bearing 
direct  and  cross-coupling  damping,  [KjJ  is  the 
bearing  axial  and  lateral  cross-coupling  stiff¬ 
ness;^^  {W  —  Wj.}  is  the  casing  vibration;  [KJ  is 
the  shaft  bow  stiffiiess;  {W  —  W^}  is  the  shaft 
residual  bow;  {F(t)}  is  the  external  and  mass 
imbalance  excitation;  and  {FQ(t)}  is  the  nonlin¬ 
ear  gear  force  through  gear-mesh  coupling. 


For  a  multiple  gear-shaft  system,  the  equa¬ 
tions  of  motion  (1)  are  repeated  for  each  shaft. 


The  motions  of  the  individual  shafts  are  coupled 
to  each  other  through  the  gear-mesh  forces 
(Fig.  1)  and  to  the  casing  through  bearing  stiff¬ 
ness  [KJ  and  damping  [Cj,]. 


The  equations  of  motion  for  the  casing  can  be 
written  as 


|Mj(w4  +  1CJ{W.  -  W}  +  (KJ{W.  -  W 
+  |C.|(W.}  +  IKjfwJ  =  (F.(t) 

(3) 

where  [MJ  is  the  mass  matrix  of  the  casing  struc¬ 
ture;  {W^}  represents  the  generalized  displace¬ 
ment  vector  of  the  casing  structure: 


(Xc) 

(Xc) 

(Yc) 

(Yes) 

(Zc) 

(Zes) 


(4) 


[CJ  is  the  casing  structure  damping  matrix;  [K^] 
is  the  casing  structure  stiffness;  and  {F^(t)}  is 
the  force  acting  on  the  casing  structure.  The 
nonlinear  gear  forces  for  the  individual  gear- 
shaft  system,  using  nonlinear  gear  stifihess^^ 
gear-tooth  friction,  are  given  as 
x-force 


®’Gxk  -  5^  Xjjji  [-Rci^ci  “  ^ck^ck 

i=l,  i  *k 

+  (Xci  -  Xck)cos  ttu  (5) 

+  (Y^  -  ou] 

[co»  «„  +  (•ipi)((i)(mii  o^)] 

y-force 
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^Gyk  -  ^tki  [“*^ci^ci  “  *^ck®ck 

i=l,i#k 

+  (Xci  -  X^)cos  ojji  (6) 

+  (Vd  -  Vd.)»ta  •»,] 

[sin  oy  +  (sign)(n)(co8  <>u)] 


where  u  is  the  rotor  critical  speed.  Similarly,  a 
set  of  orthogonality  conditions  can  be  derived  for 
the  casing  equation  of 

[MJ{WJ  +  1CJ{WJ  +  |KJ{WJ  =  0  (12) 

with  the  orthonormal  mode  [$  J  such  that 

=  II]  (IS) 

W^|CJW  =  |CJ  (1<) 


and  torsional 

®'Gtk  =  12  *^ck{*^tki  [(“^^ci^ci  “  ^ck^ck) 

i=l,i  #k 


where  is  the  casing  critical  speed.  Using 
modal  transformation^'  by  letting 


+  (Xci  Xck)co8  ttu 

+  (Y.i  -  Yi)sin  «uf 
(2) 


where  R  is  the  pitch  radius  of  the  gear  and  /i 
is  the  coefhcient  of  friction  between  the  gear 
tooth  surface  and  “SIGN”  is  the  unity  sign  func¬ 
tion  to  provide  the  sign  change  when  the  mating 
teeth  pass  the  pitch  point. 


{W} 


[^xl  {A} 
*x«l{A} 

[*y]  {B} 

[♦.]  {D} 

*t]  m 

\  / 


(16) 


Modal  Transformation  and 

The  equation  of  motion  for  the  undamped 
rotor  system  is 

(M|{W}  +  IlKJ  +  lK*jj{W)  =  0  (8) 

with  the  averi^e  bearing  support  stiffiiess  from 
the  x-y  direction^  given  as 

[Ka)  =  ^(iKtal  +  IKbyl}  (») 

M 


l^cx]  {Ac} 

(♦cx^KAc) 


{We) 


[♦cy]  {Be} 
[♦cy#l{Bc} 


The  orthogonality  condition  for  the  orthonormal 
modes  f#]  are 


(•]T(m1[*1  =  (1) 


(10) 


[*«]  {Dc} 


Pc.#l  {Dc}, 


(17) 


and 


l♦l''[K.  +  KAli»l  =  l"’l 


the  equations  of  motion  for  the  rotor  (Eq.  (1)) 
(11)  ^  transformed  as 
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[II  {Z}  +  [^,[{2}  +  [G^{Z}  +  [CJ{Z} 

-  Pb  -  Ka1<Z>  +  l«l''lCbl[».UZio}  (18) 

+  [..'|{Z}  -  [♦]■'' [Kb][«J{ZJ 
=  f»l’'{F(t)  +  FgW  +  F.(t)} 

and  for  the  casing  (Eq.  (3))  as 

|I|{Zc>  +  PJ{Z.}  +  [b-’]{ZJ 

+  [^bl^^c)  +  [1^1  {zj  -  [*J^  (19) 

WltlW  -  l»J''[Cbl[*]<z} 

=  It/iFcCt)} 


where 


The  gear  noise  rig  (Fig.  2)  was  used  to  mea¬ 
sure  the  vibration,  dynamic  load,  and  noise  of  a 
geared  transmission.  The  rig  features  a  simple 
gearbox  (Fig.  3)  containing  a  pair  of  parallel  axis 
gears  supported  by  rolling  element  bearings.  A 
150-kW  (200-hp),  variable-speed  electric  motor 
powers  the  rig  at  one  end,  and  an  eddy-current 
dynamometer  applies  power-absorbing  torque  at 
the  other  end.  The  test  gear  parameters  are 
given  in  Table  1. 

Two  sets  of  experiments  were  performed  on 
the  gearbox;  (1)  experimental  modal  analysis  and 
(2)  dynamic  vibration  measurements  during  oper¬ 
ation.  In  the  experimental  modal  analysis,  modal 
parameters,  such  as  system  naturid  frequencies 
and  their  corresponding  mode  shapes,  were  ob¬ 
tained  through  transfer  function  measurements 
by  using  a  two-channel,  dynamic  signal  analyser 
and  modal  analysis  software.  For  this  experi¬ 
ment,  116  nodes  were  selected  on  the  gearbox 


housing.  The  measured  three-dimensional  mode 
shapes  are  presented  in  Fig.  4.  The  dynamic 
vibration  measurements  consist  of  data  collected 
from  accelerometers  placed  at  three  of  the  nodes 
on  the  surface  of  the  casing.  They  were  chosen 
such  that  vibration  was  measured  in  all  three 
directions,  X,Y,  and  Z.  A  dynamic  signal  analyz¬ 
er  was  used  to  compute  the  frequency  spectra  of 
the  vibration.  The  experimental  frequency  spec¬ 
tra  are  shown  in  Figs.  5(a),  6(a),  and  7(a)  for  the 
X-,  Y-,  and  Z-directions,  respectively,  at  a  num¬ 
ber  of  different  operating  speeds. 

Discussion  and  Correlation  of  Results 

The  measured  mode  shapes,  shown  in  Fig.  4 
represent  the  major  vibration  modes  of  the  gear 
noise  rig  in  the  0-  to  3-kHz  region.  Although 
these  modes  are  only  a  small  portion  of  the  total 
modes  of  the  system,  they  represent  the  major 
portion  of  the  total  global  vibration  of  the  sys¬ 
tem.  In  order  to  produce  a  compatible  analytical 
simulation  of  the  test  apparatus,  a  similar  set  of 
modes  was  predicted  by  a  finite-element  model  of 
the  gearbox  structure.  This  model  serves  as  the 
basis  for  predicting  casing  vibrations  in  the  over¬ 
all  global  dynamic  model.  Of  the  25  modes  found 
by  the  analytical  model  in  the  0-  to  3-kHz- 
frequency  region,  the  8  dominant  modes  were 
used  to  represent  the  gearbox  dynamic  character¬ 
istics.  These  simulated  modes  are  shown  in 
Fig.  8.  The  natural  frequencies  of  the  predicted 
modes  are  within  5  percent  of  the  measured 
modes  (Table  2).  Also,  the  predicted  mode 
shapes  are  very  similar  to  the  experimental 
modes  shapes  (Fig.  4).  The  good  correlation  of 
between  the  analytical  model  and  the  measure¬ 
ments  confirms  the  accuracy  of  the  dynamic 
representation  of  the  test  gearbox  using  only  a 
limited  number  of  modes. 

For  the  dynamic  study  of  the  gearbox  vibra¬ 
tion,  it  was  found  that  during  a  slow  roll  (low- 
speed  run)  of  the  gear-rotor  assembly,  a  substan¬ 
tial  residual  bow  (or  eccentricity  in  the  sleeve 
assembly)  exists  in  the  rotor  system  as  shown  by 
its  large  orbital  motion  in  Fig.  9.  Figure  9(a) 


repr^ents  the  orbit  of  the  driver  rotor  at  the 
gear  location,  and  Fig.  9(b)  represents  that  of  the 
driven  rotor.  Note  that  the  circular  orbit  in  the 
driver  rotor  at  low  speed  represents  the  residual 
bow  deformation  of  the  rotor.  The  elliptical  orbit 
in  the  driven  rotor  is  attributed  to  a  combination 
of  the  residual  bow  effects  and  the  vertical  gear 
force  from  the  torque  of  the  driving  rotor.  In 
order  to  analytically  simulate  the  influence  of  this 
effect,  a  residual  bow  of  2  mils  (0.05  nun)  is 
incorporated  into  the  numerical  model  (Eq.  (1)). 

The  frequency  spectra  of  the  analytically  pre¬ 
dicted  casing  vibration  in  the  X-,  Y-,  and  Z- 
directions  are  presented  in  Figs.  5(b),  6(b),  and 
7(b),  respectively.  As  seen  in  Fig.  5(a),  the  ex¬ 
perimental  casing  vibration  in  the  X-direction 
shows  a  major  vibration  component  at  the  gear- 
mesh  frequency  (28  times  shaft  speed)  at  each 
rotational  speed.  A  closer  examination  of  this 
excited  frequency  component  shows  that  two 
major  vibration  peaks  occur  at  the  running 
speeds  of  1500  rpm  (at  a  tooth  pass  frequency  of 
700  Hz)  and  5500  rpm  (at  a  tooth  pass  frequency 
of  2560  Hz).  These  pedes  are  a  result  of  the 
tooth  pass  frequency  exciting  two  of  the  major 
natural  frequencies  of  the  housing,  namely  the 
658-  and  2536-Hz  modes.  The  presence  of  other 
modes  can  be  seen;  however,  the  658-  and 
2536-Hz  modes,  when  excited  by  the  gear-mesh 
frequency,  dominate  the  spectra. 

Comparing  the  predicted  vibration  spectra 
with  the  measured  spectra  reveals  that,  although 
the  actual  amplitude  values  did  not  always  agree, 
the  general  trends  of  the  spectra  were  very  simi¬ 
lar.  The  predicted  vibration  spectra  of  the  hous¬ 
ing  in  the  X-direction  is  shown  in  Fig.  5(b).  A 
comparison  of  Figs.  5(a)  and  (b)  shows  that  the 
predicted  amplitude  at  the  gear-mesh  frequency 
at  1500  rpm  is  only  3  percent  above  the  measured 
value.  The  comparison  at  5500  rpm  is  not  that 
close;  the  predicted  amplitude  is  38  percent  below 
the  measured  value.  If  trends  are  compared,  the 
predicted  spectra  show  the  same  gear-mesh, 
frequency-induced  excitation  of  the  658-  and 
2536-Hz  modes  as  that  found  in  the  measured 


spectra  at  the  running  speeds  of  1500  and 
5500  rpm,  respectively. 

Figures  6  and  7  compare  the  predicted  mea¬ 
sured  housing  vibration  spectra  in  the  Y-  and  Z- 
directions,  respectively.  The  results  of  the  com¬ 
parison  are  the  same  as  those  presented  for  the 
housing  vibration  in  the  X-direction  (Figs.  5(a) 
and  (b)).  Acutal  values  of  the  components  in  the 
spectra  were  not  always  in  good  agreement;  how¬ 
ever,  the  general  trends  for  the  predicted  and 
measured  housing  vibration  spectra  were  very 
similar.  Also,  as  seen  in  Fig.  7(b)  at  the 
1500-rpm  speed,  the  model  predicts  the  second 
and  third  harmonics  of  the  gear-mesh  frequency. 
As  shown  in  Fig.  7(a),  the  measured  vibration 
confirms  the  presence  of  these  two  harmonics  at 
the  1500-rpm  running  speed. 

Conclusions 

A  newly  developed  global  dynamic  model  was 
used  to  simulate  the  dynamics  of  a  simple 
transmission  system.  Predicted  casing  vibrations 
were  compared  to  measured  results  from  the  gear 
noise  test  rig  at  the  NASA  Lewis  Research  Cen¬ 
ter.  The  conclusions  of  thb  study  are  summa¬ 
rized  as  follows: 

1.  The  dynamics  of  the  housing  can  be  accu¬ 
rately  modeled  with  a  limited  amount  of  analyti¬ 
cally  predicted,  experimentally  verified  vibration 
modes  of  the  structure. 

2.  The  global  dynamic  model  is  capable  of 
including  in  the  analysis  the  effects  of  shaft  resid¬ 
ual  bow  or  eccentricity. 

3.  Absolute  values  of  the  housing  vibration 
predicted  by  the  global  dynamic  model  did  not 
always  agree  with  measured  values. 

4.  The  characteristics  and  trends  of  the  housing 
vibration  spectra  predicted  by  the  global  dynamic 
model  are  the  same  or  very  similar  to  those  found 
in  the  experimental  data. 
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TABLE  I.  -  TEST  GEAR  PARAMETERS 


Gear  type . Standard  involute,  full-depth  tooth 

Number  of  teeth  . 28 

Module,  mm  (diametrial  pitch  in.'^)  . 3.174(8) 

Face  width,  mm  (in.) .  6.35(0.25) 

Pressure  angle,  deg  . 20 

Theoretical  contact  ratio . 1.64 

Driver  modification  amount,  mm  (in.) .  0.023  (0.0009) 

Driven  modification  amount,  nun  (in.)  .  0.025(0.0010) 

Driver  modification  start,  deg . 24 

Driven  modification  start,  deg  . 24 

Tooth-root  radius,  mm  (in.) . 1.35(0.053) 

Gear  quality .  AGMA  cleuss  13 

Nominal  (100  percent)  torque,  N-m(in.-lb) .  71.77(635.25) 


TABLE  2.— COMPARISON  OF 
EXPERIMENTAL  MEASURED 
AND  ANALYTICAL  MODELED 
NATURAL  FREQUENICES 


Experimental, 

Hs 

Analytical, 

Hs 

Difference, 

percent 

658 

658 

0 

1049 

1006 

-4.1 

1710 

1762 

3.0 

2051 

2.6 

2276 

2336 

2.6 

2536 

1536 

0 

2722 

2752 

1.1 

2962 

3012 

1.7 

Figure  2. — Picture  of  gear  noise  rig. 


Figure  3.— Test  geartrex. 


Vibration  amplituda,  inVaac^  Vibration  ampiituda,  inysac^ 


Major  modaa  of  housing 


Gaar-mesh 

frequency 

2567-H2 

gear>niash 

frequency 


Figure  5. —  X-dtoection  experimental  and  analytical  vibration  frequency  spectra  of 
thegeaitwx 


II 


Major  modM  of  housing 


Qaar-mash 

fraquancy 


2Se7-Hz 

gaar>niash 

fraquancy 


Figura  6. —  Y>diroetlon  axparimantal  and  analytical  vibration  fraquancy  apactia  of 
thagaarbox. 
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Vibration  ampUtuda,  in7sac2  Vibration  amplituda,  In  Jaac* 


Major  modas  of  housing 


(b)  Analytical. 


Figura7.— Z*dirBCtion  axparimaiTtal  and  anaiyticai  vibration  fraquancy  spactraof 
thagaarbox. 
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Mode:  1 

Mode:  2 

Frequency  ■  658 

Frequency  >1006 

BHM 


wm 


Mode:  3 

Mode:  4 

Frequency  ■  1762 

Frequency  ■  2051 

Mode:  5 

Mode;  6 

Frequency  ■  2336 

Frequency  -  2536 

Mode;  7 

FrequeiKy  ■  2752 


Mode:  8 

Frequency  ■3012 

Figure  8.—  Qeerbox  analytical  mode  ahapee. 


Displacement,  X-direction,  in. 

(a)  Driving  rotor  orbit  (b)  Driven  rotor  orbit 
Figure  S.—Orbital  motion  of  rotor  during  slow  roll. 
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